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5 : Abstract 

o ' The charged current cross section for neutrinos with energy of a few GeV is 

1^ ■ reanalysed. In this energy range the cross section for the lowest multiplicity 

exclusive channels is an important fraction of ace and the approximation 

of describing the cross section with deep inelastic scattering formulae may 

rS \ be inaccurate. Possible consequences of our reanalysis of the cross section 

C^ ' in the interpretation of the data obtained by deep underground detectors on 

z/-induced upward going muons (both stopping and passing) are discussed. 
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It is well known that measurements of the atmospheric neutrino fluxes allow to perform 
sensitive studies on neutrino oscillations. In some cases an anomaly has been observed, 
and interpreted as positive evidence for neutrino oscillations [|T],0] . Other experiments have 
instead obtained results compatible with the no-oscillation hypothesis, and gave therefore 
exclusion plots in the neutrino oscillation parameter space (A?7i^, sin^26') 0,^]. In these 
studies it is often necessary to compare experimental results with theoretical calculations 
that depend on assumptions about the neutrino fluxes and interactions. 

In this letter we want to reanalyse the charged current (CC) neutrino and antineu- 
trino cross sections with particular attention to the energy range Ei, ^ 10 GeV, and 
discuss possible consequences for the interpretation of the measurements of the atmo- 
spheric neutrino fluxes obtained by deep underground detectors. We will argue that it is 
possible to improve on the description of the cross sections used in several recent analyses 
of neutrino-induced upward going muons [^H^ including a more careful treatment of the 
lowest multiplicity channels (quasi-elastic scattering and single pion production). The 
main effect of this more careful description of the neutrino cross sections is an increase 
of the flux of low energy upward going muons. Information on the energy spectrum of 
upward going muons can be obtained from the measurement of the ratio of event num- 
bers with muons stopping in the apparatus or passing through [Q. Since with respect to 
previous calculations 0J^ we obtain a larger rate of stopping muon events and the pre- 
dicted rate of passing events is increased by a smaller amount, the excluded region in the 
neutrino oscillation parameter space obtained from the measured ratio stopping/passing 
could be significantly reduced. The interpretation of integral measurements of the muon 
flux with a low (£'™" ~ 1 GeV) energy threshold could also be modified. 

The flux of upward going muons of energy E^ > -E™™ and direction Q can be calculated 
as 

/•oo 

^M^r, fi) = J^^^^ dE, Kii^^Eu, n) n,^(^^)_^^(i?f °; E,) , (1) 

where ^,,^(77^) (£",,, fi) is the differential flux of u^ {V^) and ni,^(j^^)^^^{E'^''^;E^) is the 
average number of muons above a threshold energy _£;™™ produced by a neutrino of energy 
E^. It is given by 
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.. «^ E.) = Na r dEo ^^^{Eo, E^) [RiEo) - R{Ef^)] . (2) 
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In equation (0), A^^ is Avogadro's number and the cross section refers to neutrino- nucleon 
CC scattering. The energy Eq of the muon at the production point is weighted by a factor 
da/dEg, the relevant cross section, and by a factor R{Eq) — R{E^^^) {R{E) is the range 
in rock of a muon of energy E) that takes into account the larger effective target available 
with increasing muon energy. Note that nj, (77 )^^t depends not only on the total CC 
cross-section but also on the shape of the muon energy spectrum. In equation (|I]) we 
are implicitly assuming that the observed muons are collinear with the parent neutrinos 
{Q^ ~ Q,y) and in eq.(|D we are neglecting fluctuations in the muon energy losses [0. 



The calculated inuon flux (eq.|l]) depends on the inclusive cross-section for muon pro- 
duction. In the literature [|-|^ this cross section has been evaluated using the deep in- 
elastic scattering formalism. The deep inelastic scattering (DIS) cross section, expressed 
in terms of the usual kinematical variables y = 1 — E^/ E^ and x = Q'^/{2mj^E^y), must 
be integrated over x and y up to the kinematical limits which, neglecting the lepton mass 
0, are < x < 1 and < y < {1 + xmj^ /{2Ey))~^ . We will denote this method of 
calculation as Method 1. 

We observe that the DIS formulae are expected to be valid only for Q^ sufficiently large, 
and that using them for calculating the cross section of low energy neutrinos implies an 
extrapolation into a region where nonperturbative effects may become important. Most of 
the parametrizations of the parton distribution functions (PDF) |^ are valid only above a 
minimum Qo of several GeV^. The parton distributions for Q^ < Ql have been considered 
to be the same as at Ql, neglecting further evolution [|,|1. More recently new sets of 
parton distributions that consider the evolution down to a lower value Ql ~ 0.3 GeV^ 



have been made available 10 



As a possible improvement we suggest to consider separately the contributions of the 
exclusive channels of lowest multiplicity, i.e. quasi-elastic scattering and single pion pro- 
duction, and describe the additional channels collectively using the DIS formulae (Method 
II). We decompose therefore the CC neutrino cross section as the sum of three contribu- 
tions: 

^1/(77) = (^QEL + O'lTT + CFdIS ■ (3) 



The quasi-elastic cross section in eq.^ is calculated following |TT]. The main uncer- 
tainty is in the axial-vector form factor, for which we follow [|12[ assuming F^iQ"^) = 



— 1.25 (1 + Q"^ /Mj^ ^, with Ma = 1-0 GeV. Inclusion of nuclear effects would decrease 
(^QEL by ~ 5%. 

To determine the energy distribution of the muon produced together with a single 
pion we make the simplifying assumption of dominance of A(1232) production. The cross 
section is normalized to the results of the more complete calculation by Fogli and Nardulli 



||13|| , for a maximum mass of the pion-nucleon system Wc = 1.4 GeV. To avoid double 



counting, the deep inelastic scattering contribution ctdis is limited to the kinematic region 



where the mass of the hadronic system in the final state isW > Wc JT^- This corresponds 
in the {x, y) plane to the condition 2 m^ E^ y{l — x) < W^ — m%. 

We observe that for neutrino energies not much larger than the nucleon mass, a large 
fraction of the phase space corresponds to W < Wc- The deep-inelastic formula in this 
region does not take into account the detailed features of the dynamics, with consequences 
both on the absolute value of the cross section and on the shape on the muon spectrum. 
Moreover, the kinematical region mj^ < W"^ < {itln + w^tt)^ is unphysical. When it is 
included in anis, as in Method I, it takes (but only roughly) into account the quasi- 
elastic contribution. 

The neutrino CC cross section and its components calculated according to equation (|]), 
using the PDF of Owens for aois, are plotted in fig.l and compared with large statistics 
data for the inclusive process at high [15| and low [16| energy, and with data on one-pion 



production [|T7| and quasi-elastic scattering [|^. The exclusive channels contribute 89% 





(12%) of the cross section for E^ = 1 (10) GeV. Our model is explicitely non-scaling, 
the ratio ay/ E^ is not a constant. In the energy range Ey ~ 1-10 GeV the ratio OyjEy 
is ~ 20% higher than the value measured at higher energy. At lower energies, closer 
to the threshold for muon production, ay/ Ey drops to zero. The comparison with data 
1T9| of the neutrino CC cross sections is rather encouraging, notwithstanding the large 
experimental errors. The situation is much less satisfactory for antineutrinos, in that our 
recipe yields results systematically larger than the sparse available data pO]. Note that 



the contribution of antineutrinos to the total flux of positive and negative muons is ~ 1/3 
of the total. 

We have calculated the flux of upward going muons according to equation (|l|) and 
its differential version using different models of the neutrino fluxes [^ and different 



parametrizations P,p!0| for the nucleon structure functions in (Tdis- In flg-2 we show the 



result of a differential muon flux calculation, in which the neutrino flux from Butkevich 
et al. I^H, and the parton distributions given by Owens have been used. The solid 



curve is obtained using Method II for the cross section (the dotted curve is the partial 
contribution of (Jdis)-, while the dashed curve is obtained with Method I (with results in 
excellent agreement with other authors [§]). The fluxes obtained with the two methods 
are essentially equal for _E^ > 10 GeV, but at lower muon energies Method II gives a result 
that is considerably larger, by ~ 12% for E^ = 3 GeV and by ~ 25% for Efj^ = 1 GeV. 
The contribution of the exclusive channels in Method II amounts to ~ 60% (32%) for 
Ef, = 1 GeV (3 GeV). 

We now discuss possible consequences of this correction to the muon flux estimate. 
The 1MB collaboration [Q has measured a ratio of rates of stopping and passing upward 
going muons {Ng/Np) = 0.16 ± 0.019. They have compared this result with a detailed 
Montecarlo calculation based on Method I, using the structure functions of EHLQ 
and the neutrino flux of Volkova, corrected at low energy with the results of Lee and 



Koh ||2T|], obtaining 0.163 for the stopping/passing ratio in the absence of oscillations. It 
has been noted that in performing the ratio stopping/passing the uncertainty due to the 
absolute normalization of the neutrino flux is greatly reduced and this has been verifled 
in independent calculations 0,0. 

To calculate a prediction for the rates of stopping and passing upward going muons in 
a speciflc detector, one needs of course a detailed knowledge of acceptances and detection 
efficiencies. As an approximation to the real experimental situation, following 0, we 
define stopping muons those in the energy interval 1.25 < E^ < 2.5 GeV and passing 
muons those with E^ > 2.5 GeV, assuming moreover that the detector acceptance and 
efficiency for muons above 1.25 GeV are approximately independent from the direction 
and energy of the particles. The fluxes $s(<l>p) are obtained integrating eq. |l|, with the 
appropriate -Emin, on the entire downward hemisphere. 

We report in Table I the results of calculations of <I>s, $p and their ratio obtained 
with different models for the neutrino fluxes and different sets of leading-order parton 
distributions. Both Method I (only DIS) and Method II - our preferred one - have been 
used. We observe that the results of Method II are larger than those of Method I. For 



those parton distributions that have a high Q^ the ratio $s/*^p is increased by ~ 10%. The 
variation is much less for parton distributions having low Q^, but it is interesting to note 
that Method II calculations using different structure functions are in better agreement 
with each other. The neutrino fluxes of Volkova and Butkevitch (and Mitsui too) predict 
very similar $s/$p, notwithstanding the differences of ~ 10% in normalization, while the 
Bartol flux is flatter and predicts a ratio somewhat smaller (by ~ 5%). 

In order to discuss the possible effects of neutrino oscillations, we define the quantity 

r = (^-/^^) (4) 

where {(^s/^p)o is the stopping/passing ratio in absence of neutrino oscillations. The ratio 
Tg depends on the neutrino oscillation parameters (Am^, sin^2^). We will only consider 
here the case of v^ <-> v^ oscillations. 

To take into account the effects of oscillations we have to make in eq. (^ the sub- 
stitution (j)y^ ^ [1 — Pv^^^ut] (pu^, and analogously for antineutrinos. The double ratio r^ 
may take the values rmin < ^s < 1- With our definitions of $s and $p, the maximum 
suppression of the stopping to passing ratio is r^i^ — 0.66, corresponding to sin^ 29 = 1 
and Am^ ~ 1.3 x 10^'^ eV^. The ratio r^ becomes unity in the limit sin^ 26* — >■ and/or 
Am^ — > 0, that correspond to the no oscillation case, and also for Am"^ ^ oo. In this limit 
(in practice for Am^ ^ 1 eV^) neutrinos of all significant energies oscillate many times, 
so that the spectrum is suppressed without distortions by a constant factor 1 — | sin^ 29, 
and the stopping/passing ratio remains unchanged. 

In fig. 3 we have drawn lines of constant r^ in the Am^,sin^26' plane. The region 
favoured at 90% c.l. by the recent Kamiokande-III analysis p[ is also shown, limited by 
the dashed line. As one can see, parameters in this region imply measurable effects in 
observations of upward going muons: in fact the region to the right of the curve r^ = 0.8 
corresponds approximately to the region excluded by the (Method I) analysis of the 1MB 
collaboration [Q. 

Let us assume that r* be the ratio of an experimental result with the theoretical no- 
oscillation value and its error Ar* includes both the experimental error and the systematic 
uncertainty in the theoretical calculation. At 90% c.l. one could exclude in the (Am^, 
sin^ 29) plane the region corresponding to values jr^ — r*| > 1.64 Ar* |^2[. Use of Method 
II for the neutrino cross section increases the theoretical prediction by ~ 10%, therefore 
r* would be lowered by the same amount, with obvious consequences on the excluded 
region. We urge therefore the experimental groups that have collected and are collecting 
data on upward going muons to reanalyze them going beyond the DIS approximation for 
the z/ cross section. 
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TABLES 

TABLE I. Fluxes (in units IQ-^^ (cm^ssr)-^) of 'stopping' (1.25 < E^ < 2.5 GeV) and 
'passing' (Ef^ > 2.5 GeV) upward going muons calculated with different models of the neutrino 
fluxes, different choices of the parton distribution functions and using Method I and Method II 
to describe the neutrino cross section. 
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FIGURES 
FIG. 1. f^ CC cross sections plotted as a function of energy. 

FIG. 2. Differential flux of upward going muons averaged in angle over one hemisphere and 
plotted as a function of muon energy. 

FIG. 3. Curves in the (A?Ti^, sin 26) plane that correspond to constant values for the 
double ratio r^. Also shown (dashed) is the 90% c.l. curve obtained by the Kamiokande-III 
combined analysis of contained and semi-contained events O] . 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411341vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411341vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9411341vl 
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Figure 2 
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Figure 3 



